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Abstract 

We consider the problem of numerical approximation for forward-backward stochastic dif- 
ferential equations with drivers of quadratic growth (qgFBSDE). To illustrate the significance 
of qgFBSDE, we discuss a problem of cross hedging of an insurance related financial deriva- 
tive using correlated assets. For the convergence of numerical approximation schemes for 
such systems of stochastic equations, path regularity of the solution processes is instrumen- 
tal. We present a method based on the truncation of the driver, and explicitly exhibit error 
estimates as functions of the truncation height. We discuss a reduction method to FBSDE 
with globally Lipschitz continuous drivers, by using the Cole-Hopf exponential transforma- 
tion. We finally illustrate our numerical approximation methods by giving simulations for 
prices and optimal hedges of simple insurance derivatives. 
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1 Introduction 

Owing to their central significance in optimization problems for instance in stochastic finance 
and insurance, backward stochastic differential equations (BSDE), one of the most efficient 
tools of stochastic control theory, have been receiving much attention in the last 15 years. A 
particularly important class, BSDE with drivers of quadratic growth, for example, arise in 
the context of utility optimization problems on incomplete markets with exponential utility 
functions, or alternatively in questions related to risk minimization for the entropic risk 
measure. BSDE provide the genuinely stochastic approach of control problems which find 
their analytical expression in the Hamilton- Jacobi-Bellman formalism. BSDE with drivers 
of this type keep being a source of intensive research. 

As Monte-Carlo methods to simulate random processes, numerical schemes for BSDE 
provide a robust method for simulating and approximating solutions of control problems. 
Much has been done in recent y e ars to cr eate schem es for BSDE with Lipschitz continuous 



drivers (see lBouchard and Touzil (|2004f ) or Elic (2006) and references therein). The numerical 



approximation of BSDE with drivers of quadratic growth (qgBSDE) or systems of forward- 
backward stochastic equations with drive r s of t his kind (qgFBSDE) turned out to be more 
complicat ed. Only recently, in dos Reis ( 2009h . one of the main obstacles was overcome. 
Following Bouchard and Touzi ( 20041) hi the setting of Lipschitz drivers, the strategy to 
prove convergence of a numerical approximation combines two ingredients: regularity of the 
trajectories of the control component of a solution pair of the BSD E in the L 2 -s ense, a 



tool first investigated in the framework of globall y Lipschitz BSDE by 
conve ni ent a priori estimate for the sol ution . See iBouchard and Touzi 



Zhand d200ll). and a 



(2004), Gobetetal 



( 2005 ). Delarue and Menozzil ( 2006 ) or Bender and Denkl (|2007fl for numerical schemes of 
BSDE with globally Lipschitz continuous dri vers, and an implementation of these ideas. 
The main difficulty treated in dos Reis ( 20091 ) consisted of establishing path regularity for 
the control component of the solution pair of the qgBSDE. For this purpose, the control 
component, known to be represented by the Malliavin trace of the other component, had 
to be thoroughly investigated in a subtle and complex study of Malliavin derivatives of 
solutions of BSDE. This study extends a thorough investigation of smoothness of systems 
of FBSDE b y methods based on Malli avin' s calculus and BMP martinga les independently 
conducted in lAnkirchner et all ( 2007b ) and Briand and Confortola ( 2008t ) . The knowledge 
of path r egular ity obtained this way is implemented in a second step of the approach in 
dos Reis (2009). The quadratic growth part of the driver is truncated to create a sequence 
of approximating BSDE with Lipschitz continuous drivers. Path regularity is exploited 
to explicitly capture the convergence rate for the solutions of the truncated BSDE as a 
function of the truncation height. The error estimate for the truncation, which is of high 
polynomial order, combines with the ones for the numerical approximation in any existent 
numerical scheme for BSDE with Lipschitz continuous drivers, to control the convergence of 
a numerical scheme for qgBSDE. It allows in particular to establish the convergence order 
for the approximation of the control component in the solution process. 

An elegant way to avoid the difficulties related to drivers of quadratic growth, and to 
fall back into the setting of globally Lipschitz ones, consists of using a coordinate transform 
well known in related PDE theory under the name "exponential Cole-Hopf transformation" . 
The transformation eliminates the quadratic growth of the driver in the control component 
at the cost of producing a transformed driver of a new BSDE which in general lacks global 
Lipschitz continuity in the other component. This difficulty can be avoided by some structure 
hypotheses on the driver. Once this is done, the transformed BSDE enjoys global Lipschitz 
continuity properties. Therefore the problem of numerical approximation can be tackled in 
the framework of transformed coordinates by schemes well known in the Lipschitz setting. 
As stated before, this again requires path regularity results in the L 2 -sense for the control 
compo nent of the so lution pair of the transformed BSDE. For globally Lipschitz continuous 
drivers IZhansd (|200ll ) provides path regularity under simple and weak additional assumptions 
such as ^-Holder continuity of the driver in the time variable. The smoothness of the 
Cole-Hopf transformation allows passing back to the original coordinates without losing 
path regularity. In summary, if one accepts the additional structural assumptions on the 
driver, the exponential transformation approach provides numerical approximation schemes 
for qgBSDE under weaker smoothness conditions for the driver. 

In this paper we aim to give a survey of these two approaches to obtain numerical results 
for qgBSDE. Doing this, we always keep an eye on one of the most important applications 
of qgBSDE, which consists of providing a genuinely probabilistic approach to utility opti- 
mization problems for exponential utility, or equivalently risk minimization problems with 
respect to the entropic risk measure, that lead to explicit descriptions of prices and hedges. 
We motivate qgBSDE by reviewing a simple exponential utility optimization problem result- 
ing from a method to determine the utility indiff erence price of an i nsuran ce rel ated asset i n 
a typical incomplete market situation, following lAnkirchner et aL ( 2007tJ ) and Freil ( 20091 ). 
The setting of the problem allows in particular the calculation of the driver of quadratic 
growth of the associated BSDE. After discussing the problem of numerical approximations, 
in this case by applying the method related to the exponential transform, we are finally able 
to illustrate our findings by giving some numerical simulations obtained with the resulting 
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scheme. 

The paper is organized as follows. In Scction[2]we fix the notation used for treating prob- 
lems about qgFBSDE and recall some basic results. Section[3]is devoted to presenting utility 
optimization problems used for pricing and hedging derivatives on non-tradable underlyings 
using correlated assets in a utility indifference approach. In section 0] we review smoothness 
results for the solutions processes of qgFBSDE, and apply them to show L 2 -regularity of 
the control component of the solution process of a qgFBSDE. In Section we discuss the 
truncation method for the quadratic terms of the driver to derive a numerical approxima- 
tion scheme for qgFBSDE. Section \6\ is reserved for a discussion of the applicability of the 
exponential transform in the qgBSDE setting. In Section[7]we return to the motivating pric- 
ing and hedging problem and use it as a platform for illustrating our results by numerical 
simulations. 



2 Preliminaries 

Fix T E K_|_ = [0,oo). We work on the canonical Wiener space (f2,.F, P) on which a d- 
dimensional Wiener process W = (W 1 ,--- , W ) restricted to the time interval [0,T] is 
defined. We denote by J- = {J-t)te[o.T] its natural filtration enlarged in the usual way by the 
P-zero sets. 

Let p > 2,m, d G N, Q be a probability measure on (f2, J 7 ). We use the symbol E® for 
the expectation with respect to Q, and omit the superscript for the canonical measure P. To 
denote the stochastic integral process of an adapted process Z with respect to the Wiener 
process on [0, T], we write Z *W = j Q Z s dW s . 

For vectors x — {x 1 , ■ ■ ■ ,x m ) in Euclidean space M m we denote \x\ — (YULii^) 2 )^ ■ I n 
our analysis the following normed vector spaces will play a role. We denote by 

• L p (R m ;Q) the space of J^-measurable random variables X : ft i— > M. m , normed by 
||A||l p = E Q [ |A| p ]p ; L°° the space of bounded random variables; 

• S p (M. m ) the space of all measurable processes (Yt)te[o,T] with values in R m normed by 

= E[(sup tg [ T j \Y t \j ]p; 5°°(R m ) the space of bounded measurable processes; 

H p (M m ,Q) the space of all progressively measurable processes (-£*)te[o,T] with values 



in W n normed by \\Z\\ W = E®[( / Q T \Z s \ 2 ds 



p/2 



BMO(F,Q) or BM02(J 7 ,Q) the space of square integrable ^-martingales $ with 
<I>o = and we set 



1$ II BMO(.F,Q) = SU P 



E Q [{$) T -($) T |J- T 



< oo, 



where the supremum is taken over all stopping times r S [0,T]. More details on this 
space can be found in Appendix 1. In case Q resp. T is clear from the context, we 
may omit the arguments Q or T and simply write BMO(Q) resp. BMO(T) etc; 

• D fc - p (K d ) and L M (M d ) the spaces of Malliavin differentiable random variables and 
processes, see Appendix 2. 

In case there is no ambiguity about m or Q, we may omit the reference to P m or Q and 
simply write S°° or H p etc. 

We investigate systems of forward diffusions coupled with backward stochastic differential 
equations with quadratic growth in the control variable (qgFBSDE for short), i.e. given 
x S W 71 , t G [O7T], and four continuous measurable functions 6, a, g and / we analyze 
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systems of the form 

In case there is no ambiguity about the initial state x of the forward system, we may and 
do suppress the superscript x and just write X, Y, Z for the solution components. For the 
coefficients of this system we make the following assumptions: 

(HO) There exists a positive constant K such that 6, er, : [0, T] x R m -> R m , 1 < i < d, are 
uniformly Lipschitz continuous with Lipschitz constant K, and 6(-,0) and <7i(-,0),l < 
i < d, are bounded by K. 

There exists a constant M G R + such that g : R m — > R is absolutely bounded by 
M, f : [0,T] x R m x I x l d -> 1 is measurable and continuous in (x,y,z) and for 
(t, x) G [0, T] x R m , y,y' eM. and 0, z' G R d we have 

\f(t,x,y,z)\<M(l + \y\ + \z\ 2 ), 
\f(t,x,y,z) - f(t,x,y',z')\ < Al[\y - y'\ + (1 + \z\ + \z'\)\z - z'\). 

The theory of SDE is well established. Since we wish to focus on the backward equation 
component of our system we emphasize that the relevant results for SDE are summarized in 
Appendix 3. 

Theorem 2.1 (Properties of qgFBSDE). Under (HO), the system j!}), (0) has a unique 
solution (X, Y", Z) G S 2 x 5°° x T-L 2 . The respective norms of Y and Z can be dominated from 
above by constants depending only on T and M as given by assumption (HO). Furthermore 

Z *W = I Z s dW s G BMO(P) and hence for allp>2 one has Z eH' p . 
Jo 

It is possible to go beyond the bounded terminal condition hypothesis by imposing the 
existence of all its exponential moments instead. In this case Z * W is no longer in BMO. 
As we shall see in Section 01 the BMO property of Z * W plays a crucial role in all of our 
smoothness results for systems of FBSDE. It combines with the inverse Holder inequality for 
the exponentials generated by BMO martingales to control moments of functional of the 
solutions of FBSDE. Smoothness of solutions is instrumental for instance in estimates for 
numerical approximations of solutions. 



x + [ b(s, X;)ds + f a(s, x;)dw s , 
Jo Jo 

g(X£) + J* f{s,X*,Y*,Z*)ds- J T Z*dW s . 



(1) 
(2) 



3 Pricing and hedging with correlated assets 



The pivotal task of mathematical finance is to provide solid foundations for the valuation of 
contingent claims. In recent years, markets have displayed an increasing need for financial 
instruments pegged to non-tradable underlyings such as temperature and energy indices or 
toxic matter emission rates. In the same manner as liquidly traded underlyings, securities on 
non-tradable underlyings are used to measure, control and manage risks, as well as to spec- 
ulate and take advantage of market imperfections. Since non-tradability produces residual 
risks which are innate and inaccessible to hedging, institutional investors look for tradable 
assets which are correlated to the non-tradable ones. In incomplete markets, one established 
pricing paradigm is the utility maximization principle. Upon choosing a risk preference, 
investors evaluate contingent claims by replicating according to an investment strategy that 
yields the most favorable utility value. Interplays and connections between the pricing of 
contingent claims on non-tradable underly i ngs and the the o ry of qgFBSDE were studied, 
among others, by lAnkirchner etaH (|2007ah . iMorlaisI (|2009f ). Ilmkeller et ail (|2009t) . and re- 
cently bv lFreil ( 2009 ). Based on this setup, we consider the problem of numerically evaluating 
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contingent claims based on non-tradable underlyings. This will be done by intervention of 
the exponential transformation of qgBSDE, to be introduced in Section [5] It allows to work 
under weaker assumptions than the numerical schemes for qgFBSDE based on the results 
reviewed in Section [4J and will allow some illustrative numerical simulations in Section [3 

The following toy market setup can be found in Section 4 of iFreil (l2009h . Assume d = 2, 



so that W = (W 1 ,W 2 ) is our basic two-dimensional Brownian motion. We use them to 
define a third Brownian motion W 3 correlated to W 1 with respect to a correlation coefficient 
p E [—1,1] according to 



Wf := 



f P AWl + f < a < T. 

Jo Jo 



Contingent claims are assumed to be tied to a one-dimensional non-tradable index that is 
subject to 

dR t = p(t, R t )dt + tr(t, R t )dWl , R = r > 0, (3) 

where fj,,a : [0, T] x 1 -> 1 are deterministic measurable and uniformly Lipschitz continuous 
functions, uniformly of (at most) linear growth in their state variable. The securities market 
is governed by a risk free bank account yielding zero interest and one correlated risky asset 
whose dynamics (with respect to the zero interest bank account numeraire) are governed by 

^ =a(s,R s )ds + f3{s,R s )dWl S = s > 0. (4) 



In compliance with lAnkirchner et al.1 (|2007afl . we assume that a,/3 : [0,T] x K — > R are 
bounded and measurable functions, and furthermore f3 2 (t,r) > e > holds uniformly for 
some fixed s > 0. Next, we set 

0(*,r):=^4, ( S ,r)e[0,T]xM, 
P{s,r) 

and note that the conditions on a and /3 imply that 9 is uniformly bounded. 

An admissible investment strategy is defined to be a real-valued, measurable predictable 
process A such that X 2 t du < oo holds P-almost surely and such that the family 



t stopping time with values in [0,T]| (5) 



is uniformly integrable. The set of all admissible investment strategies is denoted by A. In 
the following, let t € [0,T] denote a fixed time. Then the set of all admissible investment 
strategies living on the time interval [t, T] is defined analogously and we denote it by At- Let 
Vt denote the investor's initial endowment at time t, that is, vt is an J-" t -measurable bounded 
random variable. The gain of the investor at time s £ [t, T], denoted by G s , is subject to 
trading according to investing A into the risky asset, and therefore given by 

We focus on European style contingent claims, i.e. payoff profile s resuming the form F(Rt) 



where we assume, in accordance with lAnkirchner et al.l (|2007a|) , that F : K —> R is mea 



surable and bounded. Moreover the investor's risk assessment presumes that her utility 
preference is reflected by the exponential utility function, so given a nonzero constant risk 
attitude parameter 77, the investor's utility function is 

U(x) = -e- ,lx , 

The evolution of the investor's portfolio over the time interval [t, T] consists of her initial 
endowment v t , her gains (or losses) via her investment into the risky asset under an invest- 
ment strategy A and holding one share of the contingent claim F(Rt). Her objective is to 
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find an investment strategy such that her time-i utility is maximized, i.e. her maximization 
problem is given by 



V t F {v t ) := sup {E [U(v t + + F(R T ))\ F t ] : A G A t ) 

= exp{-^ t }sup {E [U(G^ + F(R T ))\F t ] : A G At) 



(0) 



For the sake of notational convenience, we write 



V F := V t F (0) = sup {E [U(G^ + F(R T ))\T t ] : A G At} ■ (7) 
Now pricing F(Rt) within the utility maximization paradigm is based on the identity 

V t °(v t ) = V t F (v t -p t ), 

where V®(v t ) den otes the time-t u t ility w ith initial endowm ent V t and with F = (see 
also Section 2 of Ankirchner et al. ( 2007al) and Section 3 of Frei ( 20091 ) ). According to 
this identity, the investor is indifferent about a portfolio with initial endowment vt without 
receiving one quantity of the contingent claim F(Rt) and a portfolio with initial endowment 
Vt—pt, now receiving one quantity of the contingent claim in addition. Hence pt is interpreted 
as the time-i indifference price of the contingent claim F(Rt). By the equality V t F (v t ) — 
exp {— rjv t } V F , it follows that 



1 , Vt 
Pt = - log— p, 



(8) 



which means that the indifference price does not depend on the initial endowment v t . Since 
the time-t indifference price ((HJ is fully characterized by V t ° an d V F , the f ocus now lies 
in the investigation of ([6]). In fact, Ankirchner et al. ( 2007a ) and Frei ( 20091 ) have already 
poin t ed ou t that ([7]) yields a characterization by means of a qgFBSDE. In accordance with 



Frei (l2009h . let u s denote by (G u )o<u<t the filtration generated by W 1 , completed by P-null 



sets. IFreil (|2009l )'s main ideas for rephrasing ((SJ) in terms of a qgBSDE are summarized in 
the following 



Lemma 3.1. The qgFBSDE 



f(u,r,z) 



F(R T ) + J f(u,R u ,Z u )du- J Z u dWl se[0,T], 
6 2 {u,r) 



2r) 



zp9(u, r) 



p )z , 



(9) 
(10) 



has a unique solution (Y, Z) G <S°° x T-L 2 such that V F = —e r,Yt holds P-almost surely. 

Proof. Si nce #(-, r) is unifor mly bounded and ^-predictable, the driver of satisfies the con- 
ditions oflKobvlanskil(l2000l) : thus © admits a unique solution (Y, Z) G 5°° x TL 2 . Moreover, 



Mania and Schweizerl (120051) have shown that Z * W 1 is both a BMO(J> and BMO(^)- 
martingale. See also lAnkirchner et al.l ( 2007a ). To prove the identity V F — —e~ vYt , we 
notice that 



-v(g$.+y t ) 



er^e 



-vYt -v(Y T -Y t ) - V (G*-Gl) 



= e -nYt e -v(YT-Y t ) e - V (Gi-G}) 



because G£ = 0. We then have 
exp{-ry(F T - Y t )} expj-^G* - G t A )} 

exp<;- // | / z„<m,;- I X u /3{u,R u )dW^ 



[X u a(u, R u ) - f(u, R u , Z u )] du 
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Denoting £f(M) = £(M) s /£(M) t for t < s < T where £(M) S is the stochastic exponential 
of a given semi-martingale M, we introduce 

K u ■= \ (v {pZ u + [3(u, R u )\ u ) - 9 u y, t<u<T. 
Then a simple calculation yields 

exp { - r,(Y T - Y t ) } exp { - V (G* - G$) } 

= £j (^J -nZdW 1 - J rjX/3(-, R)dW 3 ^ exp j jf if„du 1 . 



Since A/3(-, i?) * is a BMO-martingale, we can condition with respect to the cr-algebra T t 
and get 



E 



= e-^e£ K " du >e-^. (11) 



By ([5]) and a localization argument, this inequality holds for every A £ At, and therefore we 
have V t F < —e~ 7,Yt . To prove equality, note that the inequality (ITT]) becomes an equality for 
^ u = ~ fs {u,R ) ^ u ■qplu^R \ ' ^ ms m conjunction with the observation that 

exp j-^A^j = exp = exp [~ v (g* - G\) } 

= £j Q -nZdW 1 - J rf\p{; R)dW^j x exp {- v (Y T - F t )} 

is the product of a bounded process and true ^-martingale yields that condition ([5]) is 
satisfied. Hence A £ At and we have shown V t = —er"^ 1 . □ 

The proof of the previous Lemma 13.11 yields the following 
Corollary 3.2. The investment strategy 

P{s,R s ) i]/3{s,R s ) 
where Z is the control component of the solution to ®, belongs to At and satisfies 



E 



U(v t + G Xt + F{R T )\F t = sup {E [U(v t + G* + F(R T )\ Ft] ■ A 6 At} = V t F (v t ). 



One application is given in the following example 
Exam ple 3.3. [Put option on kerosene, compare with Example 1.2 from \Ankirchner et ~al. 



i < £007d) J Facing recent considerable declines in world oil prices, companies producing kerosene 
wish to partially cover their risk of such a depreciation. European put options are an estab- 
lished financial instrument to comply with this demand of risk covering. Since kerosene is 
not traded in a liquid market, derivative contracts on this underlying must be arranged on 
an over-the-counter basis. Knowing that the price of heating oil is highly correlated with the 
price of kerosene, the pricing and hedging of a European put option on kerosene can be done 
by a dynamic investment in (the liquid market of) heating oil. A numerical treatment of this 
pricing problem will be displayed in Section^ 



4 Smoothness and path regularity results 

The principal aim of this paper is to survey some recent results on the numerical approxi- 
mation of prices and hedging strategies of financial derivatives such as the liability F(Rt) 
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in the setting of the previous section. As we saw, this leads us directly to qgFBSDE. In 
the subsequent sections we shall discuss an approach based on a truncation of the driver's 
quadratic part in the control variable. It will be crucial to give an estimate for the error 
committed by truncating. Our error estimate will be based on smoothness results for the 
control component Z x of solutions of the BSDE part of our system. Smoothness is under- 
stood both in the sense of regular sensitivity to initial states x of the forward component, as 
well as in the sense of the stochastic calculus of variations. Since the control component of 
the solution of a BSDE is related to the Malliavin trace of the other component, we will be 
led to look at variational derivatives of the first order. 

Our first result concerns the smoothness of the map [0,T] x M. m 9 (t,x) \-t Z x , espe- 
cially its differentiability in x. The second result refers to the variational differentiability of 
(Y x , Z x ) in the sense of Malliavin's calculus. We shall work under the following hypothesis, 
where we denote the gradient by the common symbol V, and by V M if we wish to emphasize 
the variable u with respect to which the derivative is taken. 

(HI) Assume that (HO) holds. For any < t < T the functions b(t, •), <Ji(t, •), 1 < i < d, 
are continuously diffcrcntiable with bounded derivatives in the spatial variable. There 
exists a positive constant c such that 

y T a(t,x)a T (t,x)y>c\y\ 2 , x, y € R m , t € [0, T]. (13) 

/ is continuously partially differentiable in [x, y, z) and there exists M £ M+ such that 
for (t, x, y, z) G [0, T] x K m xRx R d 

\W x f(t,x,y,z)\<M(l + \y\ + \z\ 2 ), 
\V y f(t,x,y,z)\<M, 
\V z f(t,x,y,z)\ <M(l + \z\). 
g : R m — > M. is a continuously differentiable function satisfying \Vg\ < M . 



Smoothness results 



The fol lowin g differentiability results are ex tensions of Theorems proved in lAnkirchner et al 



(|2007bT ) and iBriand and Confortolal ([20081 ). For further de tails, comments and complete 
proofs we refer to the mentioned works or to ldos Reisl f|2009h . 

Theorem 4.1 (Classical differentiability). Suppose that (HI) holds. Then for all p > 2 the 
solution process Q x = (X X ,Y X , Z x ) of the qgFBSDE £7J), $j) with initial vector x 6 M m for 
the forward component belongs to S p x S p x H p . The application W n 3 x h-> (X x , Y x , Z x ) e 
5 p (]R m ) x S p (R) x H p (M. d ) is differentiable. The derivatives of x i-> X x satisfy $Tj) while 
the derivatives of the map x i— > (Y x , Z x ) satisfy the linear BSDE 

VY t x = Vg(X x )VX x - J VZ*dW s + J (Vf(s,Q x ),Ve x )ds, te[0,T\. 

Theorem 4.2 (Malliavin differentiability). Suppose that (HI) holds. Then the solution 
process (X, Y, Z) of FBSDE flj), fH) has the following properties. For x S R m , 

• X x satisfies (MS) and for anyO<t<T,x£ W n we have (Y x , Z x ) e L x , 2 x (Li :2 ) d . X x 
fulfills the statement of Theorem \7.4\ and a version of (D U Y X , D u Z x )o< u ,t<T satisfies 



D U Y X = 0, 



D U Z X = 0, 



D U Y X = Wg(X x )D u X x 



t <u<T, 



(vf(s,e x s ),D u e x )ds 



D U Z X S &W S , te[u,T}. 



Moreover, (D t Y x )o<t<T defined by the above equation is a version of (Z x )o<t<T- 
The following representation holds for any < u < t < T and x £ K m 

D U Y X = V x Y t x (V x X:)- l a(u, XZ), a.s., 
Z t - V X Y X {V x X x y 1 a{s,X x ), a.s.. 
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Regularity and bounds for the solution process 

A careful analysis of DY in both its variables under the smoothness assumptions on the 
coefficients of our system formulated earlier reveals the following continuity properties for 
the control process Z. 

Theorem 4.3 (Time continuity and bounds). Assume (HI). Then the control process Z 
of the qgFBSDE {2]j-0) has a continuous version on [0,T]. Furthermore for all p > 2 it 
satisfies 



\Z\\ S p < oo. 



(14) 



Theorem 4.4 (Regularity). Under (HI) the solution process (X, Y, Z) of the qgFBSDE (J\), 
(0) satisfies for all p > 2 

i) there exists a constant C p > such that for < s < t < T we have 

E[ sup \Y U -Y S \ P ] < C p \t~s\%; 

s<u<t 

ii) there exists a constant C p > such that for any partition ir — {to, ■ ■ -t^} with = 
to < • • • < ijv — T of [0, T] with mesh size \tt\ 



\Z t -z t fdt 



< C p \tt\ 



Now let h — T/N, tt n = {tj = ih : i = 0, • • • , N} be an equidistant partition of [0,T] 
with N + 1 points and constant mesh size h. Let Z be the control component in the solution 
of the qgFBSDE (fTJ), @ under (HI) and define the family of random variables 



zi 



Z s ds J*, 



U 6 ^ \ {t N }. 



(15) 



For < i < N — 1 the random variable Z*. is the best -measurable approximation of Z 
mH 2 ([U,t l+1 }), i.e. 



E 



U+i 



\z„ - z 



w |2 ds 



= infE 

A 



U+i 



\Z„ - K\ 2 ds 



where A is allowed to vary in the space of all square integrable J 7 ^ -measurable random 

— N — N 

variables. By constant interpolation we define — Z^. for t 6 ti+i [, < i < N — 1. It 
is easy to see that (Z^ N ) te [ ^ converges to (Z t ) t6 [ .T] in "H 2 [0,T] as h vanishes. Since Z is 
adapted there exists a family of adapted processes Z™ indexed by our equidistant partitions 

N N o 

such that Zf — Z ti for t £ [ti, U + i) and that Z* converges to Z in H as h tends to zero. 
Since Z v is the best "H 2 -approximation of Z, we obtain 



Z* 



\u- 



< \\Z - z n ' 



— !> 0, as h 



0. 



The following Corollary of Theorem 22] extends Theorem 3.4.3 in Zhane ( 200lh (see Theorem 
173]) to the setting of qgFBSDE. 

Corollary 4.5 (L 2 -regularity of Z). Under (HI) and for the sequence of equidistant parti- 
tions (7r A ')jveN of [0,T] with mesh size h — j£ , we have 



max 

0<i<N- 



< sup E 

1 te[ti,t i+ i) 



i=0 Jti 



- — jv ,o 

K l 2 d s 



< Ch, 



where C is a positive constant independent of N . 
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Remark 4.6. The above corollary still holds if (HI) is weakened. More precisely, the corol- 
lary's statement remains valid if one replaces in (HI) the sentence 

"g : R m — 5- K is a continuously differentiable function satisfying |Vg| < M." 

by 

"g : R m —¥ R is uniformly Lipschitz continuous in all its variables. " 
The proof requires a regularization argument. 

5 A truncation procedure 

To the best of our knowledge so far none of the usual discretization schemes for FBSDE 
has been shown to converge in the case of systems of FBSDE considered in this paper, the 
driver of which is of quadratic growth in the control variable. The regularity results derived 
in the preceding section have the potential to play a crucial role in numerical approximation 
schemes for qgFBSDE. We shall now give arguments to substantiate this claim. In fact, the 
regularity of the control component of the solution processes of our BSDE will lead to precise 
estimates for the error committed in truncating the quadratic growth part of the driver. We 
will next explain how this truncation is done in our setting. 

We start by introducing a sequence of real valued functions (hn)n&\ that truncate the 
identity on the real line. For n£N the map h n is continuously differentiable and satisfies 

• h n — > id locally uniformly \h n \ < |id| and \h n \ < n + 1; moreover 

( (n + l) ,x>n + 2, 
h n (x) = < x ,\x\<n, 

[ -(n+l) ,x < -(n + 2); 

• the derivative of h n is absolutely bounded by 1 and converges to 1 locally uniformly. 

We remark that such a sequence of functions exists. The above requirements are for instance 
consistent with 

h n (x) - 



- n 2 + 2nx - x(x - 4))/4 , x G [n, n + 2] , 



in 2 



2nx + x(x + 4))/4 , x G [-(n + 2), -n]. 



We then define h n : R d —5- R d by z h n (z) = (h n (zi), • • ■ , h n (zd)), n e N. The sequence 
(^n)neN is chosen to be continuously differentiable because the properties stated in Theorem 
14.31 need to hold for the solution processes of the family of FBSDE that the truncation 
sequence generates by modifying the driver according to the following definition. 

Recalling the driver / of BSDE (J2J, for n £ N we define f n (t, x, y, z) := f(t ,x,y, ti n 
(t, x, y, z) e [0, T] x R m x R x R d . With this driver and flU we obtain a family of truncated 
BSDE by 

Y t n =g(X T ) + £ f n (s,X s ,Y s n ,Z?)ds- £ Z«dW 8 , te[0,T],n&K (16) 

The following Theorem proves that the truncation error leads to a polynomial deviation 
of the corresponding solution processes in their natural norms, formulated for polynomial 
order 12. 

Theorem 5.1. Assume that (HI) is satisfied. Fix ngN and let X be the solution of $7\). 
Let (Y, Z) and (Y n , Z™)„ s n be the solution pairs of (0) and &16]) respectively. Then for all 
p > 2 there exists a positive constant C p such that for all n € N 

cT 



E 



sup \Y t n -Y t \p] +e\( [ \Z:-Z s \ 2 d 
te[o,T] LV Jo 



<C J- 



The proof of Theorem 15.11 roughly involves estimating the probability that Z n exceeds 
the threshold n as a function of n G N through Markov's inequality. The application of 
Markov's inequality is possible thanks to (Tl4]) . 
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6 The exponential transformation method 



In the preceding sections we exhibited the significance of path regularity for the solution of 
systems of qgFBSDE, in particular the control component, for their numerical approximation. 
In this section we shall discuss an alternative route to path regularity of solutions in a 
particular situation that allows for weaker conditions than in the preceding sections. We 
will use the exponential transform known in PDE theory as the Cole-Hopf transformation. 
This mapping takes the exponential of the component Y of a solution pair as the new first 
component of a solution pair of a modified BSDE. It makes a quadratic term in the control 
variable of the form z > j\z\ 2 vanish in the driver of the new system. The price one has to 
pay for this approach is a possibly missing global Lipschitz condition in the variable y for the 
modified driver. It is therefore not clear if the new BSDE is amenable to the usual numerical 
discretization techniques. We give sufficient conditions for the transformed driver to satisfy 
a global Lipschitz condition. In this simpler setting our techniques allow an easier access to 
smoothness results for the solutions of the transformed BSDE. The Cole-Hopf transformation 
being one-to-one, it is clear that regularity results carry over to the original qgFBSDE. 

Under (HO), we consider the transformation P = e lY and Q = jPZ. It transforms our 
qgBSDE © with driver / into the new BSDE 

P t =e?^+ [ lPs f[s,xJ^,^y\^]d S - J\ S AW S , te[0,T}. (17) 

Combining ([TTf with SDE (fTJ), we see that for any p > 2 a unique solution (X,P,Q) G 
S p x S°° x H p of (fTJ and (IT7|) exists. The properties of this triple follow from the properties 
of the solution (X,Y,Z) of the original qgFBSDE (JTJ) and ©. For clarity, we remark that 
since Y is bounded, P is also bounded and bounded away from 0. The latter property allows 
us to deduce from the BMO martingale property of Z * W the BMO martingale property of 
Q * W. For the rest of this section we denote by K. a compact subset of (S, +oo) for some 
constant 5 € K+ in which P takes its values. 

The form of the driver in ([T7|) indicates that after transforming drivers of the form of the 
following hypothesis, we have good chances to deal with a Lipschitz continuous one. 

(HO*) Assume that (HO) holds. For 7 6 K let / : [0, T] x M m x M x R d ->■ M be of the 

form 

/(*, x, y, z) = l(t, x, y) + a(t, z) + ^|z| 2 , 

where I and a are measurable, I is uniformly Lipschitz continuous in x and y, a is 
uniformly Lipschitz continuous and homogeneous in z, i.e. for c£K, (s, z) g [0, T] x R d 
we have a(s, cz) = ca(s, 2); I and a continuous in t. 

Assumption (HO*) allows us to simplify the BSDE obtained from the exponential transfor- 
mation to 



Pt = e i9{x T ) + f F ( SjXs ,P s ,Q s )ds- [ Q s dW s , t£[Q,T], 
Jt Jt 

where the driver is defined by 

F : [0, T] x JR m x K x R d -> K, 

(s,x,p,q) i-> ypl(s,x, ^if.) + lpa (s, —). (19) 
v 7 v 7p' 

Thanks to the homogeneity assumption on a our driver simplifies further. Indeed, we have 
for (s,x,p,q) e [0,T] x E m x K x R d 

F(s,x,p,q) =^pl(s,x, ] ^-) + a(s,q). (20) 

The terminal condition of the transformed BSDE still keeps the properties it had in 
the original setting. Indeed, boundedness of g is inherited by exp(7<?). Furthermore, if g is 



(18) 
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uniformly Lipschitz, then clearly by boundedness of g, the function e 19 is uniformly Lipschitz 
as well. 

Let us next discuss the properties of the driver ([TT)]) in the transformed BSDE. We recall 
that since I and a are Lipschitz continuous, there is a constant C > such that for all 
(s, x,p, q) G [0, T] x R m x K, x W l we have 

\F(s,x,p,q)\ < \jpl(s,x,-^-) +a(s,q)\ 

< C\p\(l + \x\ + | logp | + \q\) < C(l + \x\ + \ P \ + \q\) . 

This means that F is of linear growth in x,p and q. 

To verify Lipschitz continuity properties of F in its variables x,p and q, by (|20p and the 
Lipschitz continuity assumptions on a, it remains to verify that 

t \ , u lo % p \ 
(x,p) i — ^ 7pt(s, j 

7 

is Lipschitz continuous in a; and p, with a Lipschitz constant independent of s G [0,T]. As 
for a;, this is an immediate consequence of the Lipschitz continuity of I in x. For p we have 
to recall that p is restricted to a compact set K, C M+ not containing 0, to be able to appeal 
to the Lipschitz continuity of I in y. This shows that F is globally Lipschitz continuous in 
its variables x,p and q. 

We may summarize these observations in the following Theorem. 

Theorem 6.1. Let f : [0,T] x R m x M x M d — >• K be a measurable function, continuous 
onR m x M x R d , and satisfying (HO*). Then F as defined by U9\) is a uniformly Lipschitz 
continuous function in the spatial variables. 

Theorem 16.11 opens another route to tackle convergence of numerical schemes via path 
regularity of the control component of a solution pair of a qgFBSDE system whose driver 
satisfies (HO*). Look at the new BSDE after applying the Cole-Hopf transform. Since it 
possesses a Lipschitz continuous driver, path regularity for the control component Q of the 
transformed BSDE will follow from Zhang's path regularity result stated in (|7.5I) provided 
the driver is ^-Holder continuous in time. Of course, by the smoothness of the Cole-Hopf 
transform, the control component Z of the original BSDE will inherit path regularity from 
Q. This way we circumvent the more stringent assumption (HI) which was made in section 
IH 

In what follows the triples (X, Y, Z) and (X, P, Q) will always refer to the solution of 
qgFBSDE ©, © and FBSDE (TTJ), {TSJ respectively. 

Theorem 6.2. Let (HO*) hold. Assume that 

[0, T]xM m xKxR^ (s, x,p, q) ^ F{s, x,p, q) E R, 

the driver of BSDE U8\) . is uniformly Lipschitz in x,p and q and is ^-Holder continuous 
in s. Suppose further that the map g : Mr M, as indicated in (HO), is globally Lipschitz 
continuous with Lipschitz constant K . Let {X, Y, Z) be the solution of qgFBSDE {Ip, (0), and 
e > be given. There exists a positive constant C such that for any partition iz = {to , • • • ,i./v} 
with = to,T = t]y,to < • • • < tjy °f the interval [0,T], with mesh size \ir\ we have 

max { sup Ehy t -y t .| 2 j} <C\tt\ and Vlf \Z S - Z£ | 2 ds] < C^] 1 ^ . 

0<i<N-i I te[u,t i+1 ) L J J T~i L Jt, J 



i=0 



Moreover, if the functions b and a are continuously differentiable in x G R m then t h- > Z t is 
a.s. continuous in [0,T]. 

Proof. Throughout this proof C will always denote a positive constant the value of which 
may change from line to line. Let (X,P,Q) be the solution of (HJ and (fT5)l. where P takes 
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its values in K. and Q * W is a BMO martingale. Applying Theorem 17.51 yields a positive 
constant C such that for any partition n = {to, ■ ■ ■ , t^} of [0, T] with mesh size \ir\ 



max 

0<i<N- 



,{ 



sup E 

t£[ti,t i+1 ) 



\p t -p ti \ 



Since P takes its values in the compact set K, C 
C such that for any < i < iV — 1, i G [tj, 



< C|tt|. 

not containing there exists a constant 



\Y t -Y u \ = C\ log P t - log P ti | < C\P t - P u | 
Using the two above inequalities we have 

| 2 H<C max ( sup E 
1 J J o<i<w-i ItefoA+o 



max 

0<i<JV- 



{ sup E 

1 L te[u,t i+1 ) 



\Y t - Yt 



•tj 2 ]}<cw. 



This proves the first inequality. For the second one, note that by definition for < i < 

N-i,te[ti,U+i) 



< \z t 



< 



;{'f-f' + 'f-ftH{i< 

<c{|Q*| \P t -Pu\ + \Qt~Qu\}- 
We therefore have for < i < N — 1 



1 



E 



|Z S -ZZ\ 2 ds 



< E 



Z,. | 2 ds 



< 2C\E 



sup |P t -P ti | 

*e[t»,ti+i) 



\ 2 ds 



■ E 



|Q t -Q tl | 2 ds]}. 



Since Q G H p for all p > 2, for any two real numbers a, (3 G (1 , oo) satisfying 1 / a+ 1//? = 1 we 
may continue using Holder's inequality on the right hand side of the inequality just obtained, 
and then Theorem 17.51 to the term containing P. This yields the following inequality valid 
for any < i < N — 1 with a constant C not depending on % 



E 



U+i 



\Z„-Zl\ 2 ds 



< C 



sup \P t 



Pt. 



7 E[ 



\ 2 ds 



h} 



<C\E[ sup \P t -P u \ 2 Y +|tt|} <c\\nfi + \n\\. 
1 te[u,u+i) ' y 3 



Now choose a — tzt, to complete the claimed estimate. 

To prove that Z admits a.s. a continuous version, it is enough to remark th at the 
Theorem's assumptions imply the conditions of Corollary 5.6 in lMa and Zhang (|2002). The 
referred result yields that Q is a.s. continuous on [0, T]. Since P is continuous and bounded 
away from zero we conclude from the equation 7PZ = Q that Z is a.s. continuous as 
well. □ 



7 Back to the pricing problem 

We now come back to the numerical valuation of the put option on kerosene as depicted in 
example 13.31 Notations in the following are adopted from Section [3] Assume that the put 
option expires at T — 1 . Let R and S denote the dynamics for the financial value of kerosene 
and heating oil respectively. In particular we assume both dynamics to be lognormally 
distributed according to 

dR t = fi(t,R t )dt + a(t,R t )dWl = 0.12 R t dt + 0A1 R t dWl , 

A Q 

—1 = a (t, R t )dt + 0(t, R t )dW? = 0.1 dt + 0.35 dW t 3 , 
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and we assume the spot price for heating oil to be sq — 173 money units (e.g. US Dollar, 
Euro), see also equations Q and (0]). Risk aversion is set at the level of -q = 0.3. Figure 
H] displays sample paths of the kerosene price with a spot price of r = 170 and heating 
oil price at different correlation levels using the explicit solution formula for the geometric 
Brownian motion. We see that the higher the correlation, the better the approximation of 
the kerosene by heating oil becomes. We have seen that the valuation of the put option via 



- # - heating oil, rho-0.99 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 1: Price paths of the nontradable asset kerosene and the correlated asset heating oil at 
different correlation levels. The spot of kerosene was set to vq = 170. 

utility maximization yields the pricing formula ([5]) which in conjunction with Lemma 13.11 
becomes the difference of two solutions of a qgBSDE with the generator (flUl) 

Pt=Y t F -Y t °, 0<i<T, 

where F(x) = (K — x) + for some strike K > 0. For the numerical simulation of the qgFBSDE 
Y F and Y°, we apply the exponential transformatio n to both BSDE (see Section[B]) and then 
employ the algorithm by iBender and Denk ( 2007 ) with N — 100 equidistant time points, 
70000 paths and a regression basis consisting of five monomials and the payoff function of the 
put option. The Picard iteration stops as soon as the difference of two subsequent time zero 
values is less than 10~ 5 . Simulations reveal that 12 to 13 iterations are needed for solving 



one exponentially transformed qgFBSDE. Figures 2(a) and 2(b) depict the time zero price 
Po of the put option at different strike and kerosene spot levels. The lower the correlation, 
the lower the price becomes. This is clear because lower correlations between heating oil 
and kerosene lead to higher non-hedgeable residual risk which diminishes the risk covering 
effect of the contingent claim and thus also its value. Figures 3(a) and |3(b)] depict sample 
paths of the dynamics for the price p t and the optimal investment strategy 7r t for an at the 
money put with strike K = 180 and kerosene spot r$ — 170. The plots depict price and 
monetary investment for every fourth time point of the discretization. The price process 
and the dynamics of the optimal investment strategy are intertwined: high fluctuations of 
the price process result in high fluctuations of the investment strategy and vice versa. In 
general we observe that replication on high correlation levels tends to entail greater market 
activity because kerosene price risks can then be well hedged by market transactions that 
move closely along the dynamics of heating oil. In contrast, replication on lower correlation 
levels leads to a higher amount of residual risk which is inaccessible for hedging and thus 
lower market activity is needed. 
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(a) Put option price in terms varying strikes at a fixed (b) Put option price in terms of varying kerosene spots 
kerosene spot ro = 170. at a fixed strike K — 200. 

Figure 2: Values of the put option in terms of kerosene spot and strike for varying correlations. 
High correlations lead to high the prices for the contingent claim. 
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- - rho=0.97 
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0.1 0.2 0.3 0.4 0.5 0.6 
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-©-rho=0.97 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
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(a) Dynamics of the price process pt for strike K 
180. 



(b) Dynamics of the optimal investment strategy 7Tt 
for strike K = 180. 



Figure 3: Paths of the price pt and the optimal investment strategy irt for varying correlation 
levels. In general high correlations entail greater market activity. 



Appendix 1 — Some results on BMO martingales 



BMO martingales play a key role for a priori estimates needed in our s ensitivity analysi s of 
solutions of BSDE. For details about this theory we refer the reader to iKazamakil (|199J). 

Let $ be a BMO(J-, Q) martingale with $o = 0. $ being square integrable, the 
martingale representation Theorem yields a square integrable process </> such that <I>t = 
J*4> a dW a ,t G [0,T]. Hence the BMO(F,Q) norm of $ can be alternatively expressed as 



sup K- 

t J~ — stopping time in [0,T] 



1 As\Tt 



< 00. 



Lemma 7.1 (Properties of BMO martingales). Let $ be a BMO martingale. Then we have: 

1) The stochastic exponential £(<&) is uniformly integrable. 

2) There exists a number r > 1 such that £(<&t) £ L r . This property follows from the 
Reverse Holder inequality. The maximal r with this property can be expressed explicitly 
in terms of the BMO norm of $. 
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3) If <f> = J Q 4> s ds has BMO norm C , then for p > 1 £/ie following estimate holds 

E[(^ T |^| 2 d S ) P ]<2p!(4C 2 f. 
iTence BMO C H p /or allp>\. 



Appendix 2 — Basics of Malliavin's calculus 

We briefly introduce the main notation of the stochastic calcu lus of variation s also known as 
Malliavin's calculus. For more details, we refer the reader to iNualarti ((2006). Let S be the 
space of random variables of the form 

£ = f([J^ hlMw^^, ■■■,(£ h d /dWf )i<*<„)) , 

where F G C b °°(M Tlxd ), h 1 ,-- - , h n G L 2 ([0, T]; R d ), n G N. To simplify notation, assume that 
all/i 3 are written as row vectors. For £ G 5, we define D = (D 1 , • • • ,D d ) : S — > L 2 (Qx [0, T]) d 

by 



1=1 1 J0 



h n t m t )h^\ < 9 < T, 1 < i < d, 



and for k G N its fc-fold iteration by 

For fc G N, p > 1 let D fe ' p be the closure of iS with respect to the norm 

iieii^Epieii^+Eiii^^niW 

i=i 



is a closed linear operator on the space H) k > p . Observe that if £ G D 1 ' 2 is .Fj -measurable 
then D e £ = for 9 G (t,T]. Further denote EP fe >°° = n p> iD fc ' p . 

We also need Malliavin's calculus for smooth stochastic processes with values in R m . For 
k G N,p > 1, denote by Lfc iP (R m ) the set of R m -valued progressively measurable processes 
it = (u x , ■ ■ • , u m ) on [0, T] x ft such that 

i) For Lebesgue-a.a. t G [0,T], •) G (D fc 'P) m ; 

ii) [0,T]xfi 9 (t,u) h-> D {kS> u(t,uj) G (L 2 ([0, T] 1+fe )) <ixn admits a progressively measurable 
version; 

Note that Jensen's inequality gives for all p > 2 

/.T 



E 



o Jo 



\D u X t \ 2 dudt 



\D u Xf np du. 



Appendix 3 — Some results on SDE 

We recall results on SDE known from the literature that are relevant for this work. We state 
our assumptions in the multidimensional setting. However, for ease of notation we present 
some formulas in the one dimensional case. 
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Theorem 7.2 (Moment estimates for SDE). Assume that (HO) holds. Then j7J) has a 
unique solution X G S 2 and the following moment estimates hold: for any p > 2 there exists 
a constant C > 0, depending only on T, K and p such that for any x £ W 71 , s, t £ [0, T] 



E[ sup |Xt| p ] < CE 

0<t<T 



(|6(t,0)| p + k(t,0)|f)di 



E[ sup \X U -X s \ p ] < CE sup {|6(i,0)| p + |(j(t,0)| p } |i-s| p/2 . 

Furthermore, given two different initial conditions x, X £ W 71 , we have 

E[ sup pff-JCf'H <C|a;-x'| p . 

Theorem 7.3 (Classical differentiability). Assume (HI) holds. Then the solution process X 
of (OP as a function of the initial condition x £ R m is differentiate and satisfies for t £ [0, T] 



VX t = J ro + / V6(X s )VX s ds + f Va(X s )VX s dW s 
Jo Jo 



(21) 



where I m denotes the m x m unit matrix. Moreover, VX t as an m x m-matrix is invertible 
for any t £ [0, T]. Its inverse (VX 4 ) _1 satisfies an SDE and for any p > 2 there are positive 



constants C p and c p such that 



IVXI 



|(vx)- 



5f 



<c„ 



and 



E 



sup \(VX u )-(VX s )\p + sup KVX,,)- 1 - (V^)- 1 !" <c p \t-s\r/ 2 . 

s<u<t s<u<t ^ 



Theorem 7.4 (Malliavin Differentiability). Under (HI), X £ L^2 <™d its Malliavin deriva- 
tive admits a version (u, t) h-» D u X t satisfying for < u < t < T the SDE 



D u X t = a(X u )+ f Vb(X s )D u X s ds+ f Wa(X s )D u X s dW s 



(22) 



Moreover, for any p > 2 there is a constant C p > such that for x £ R m and < v < u < 
t<s<T 

\\D u X\\ p SP <C p (l + \xf), 
E[ \D u X t - D U X S \P] <C P (1 + \x\P)\t -s\i, 
\\D U X - D V X\\% < C p (l + \x\ p )\u -v\*. 

By Theorem \ 7. 3\ we have the representation 

D u X t = VX t (VX u )- 1 cr(X 1 ,)i [ o,t,](t), for all u,t £ [0,T]. 



Appendix 4 — Path regularity for Lipschitz FBSDE 



We state a version of the L 2 -regularity result for FBSDE satisfying a global Lipschitz condi- 
tion. The result which was seen t o be closely re lated to the convergence of numerical schemes 
for systems of FBSDE is due to IZhand (|200lh . For our FBSDE system © we assume 
that 6, cr, /, g are deterministic measurable functions that are Lipschitz continuous with re- 
spect to the spatial variables and ^-Hol der continuous wi t h resp ect to time. Furthermore we 
assume that a satisfies (fT5|) . Then from lEl Karoui et~all (Il997h one easily obtains existence 
and uniqueness of a solution triple (X, Y, Z) of FBSDE (TTJ), © belonging to S 2 x S 2 x H 2 . 
For a partition 7r of [0, T] define the process as in ([Poll . Then the following result holds. 
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Theorem 7.5 (Path regularity result of lZhand (|200lh ). Let (X, Y, Z) £ S 2 x S 2 x H 2 be the 
solution of FBSDE (QP 7 (0) in the setting described above. Then there exists C € R+ such 
that for any partition tt = {to, • • • , ijv} of the time interval [0, T] with mesh size \ir\ we have 



max < sup E 

0<i<JV-l l te[tj ,t i+1 ) 

N-l r U+i 



\Yt - Y, 



i=0 



|Z s -Z£| 2 d S 



w-i 



i=0 



ti+1 



<C\ir\. 
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